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Recently exact numerical results from the evolution equation for parton multiplicities in QCD jets have been 
obtained. A comparison with various approximate results is presented. A good description is obtained not only 
of the jet multiplicities measured at LEP-1 but also of the hadron multiplicities for cms energies above 1.6 GeV 
in e^e~annihilation. The solution suggests that a final state hadron can be represented by a jet in the limit of 
small (nonperturbative) fcx cut-off Qo. In this description using as adjustable parameters only the QCD scale A 
and the cut-off Qo, the coupling as can be seen to rise towards large values above unity at low energies. 



1. INTRODUCTION 

The phenomena of multiparticle production in 
hard processes require for their description within 
QCD both perturbative and non-perturbative el- 
ements. Whereas the perturbative part - at least 
in principle - is rather well understood, there are 
different approaches to deal with the hadroniza- 
tion at large distances and also the characteristic 
scale for the transition between both regimes is 
not uniquely defined. Recently, we studied this 
transition region more quantitatively in case of a 
simple observable, the multiplicity of jets in an 
event 

At a large resolution scale Qc only a few jets 
are resolved and their rate can be calculated in 
fixed order perturbation theory; with increasing 
resolution more and more jets are resolved and 
in this case results can be obtained by a resum- 
mation of the leading and next-to-leading loga- 
rithmic terms of the perturbation theory [Q. In 
these calculations the multiplicity is obtained in 
absolute normalization and the only free param- 
eter is the QCD scale A. Increasing the resolu- 
tion further the jet multiplicity coincides finally 
with the hadron multiplicity. The jet multiplicity 
is defined for a given resolution cut-off parame- 
ter Qc using the well known Durham/A: algo- 
rithm which counts some clusters of particles as 
separate if their relative transverse momenta are 
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larger than Qc- In e'''e~one also introduces the 
normalized scale yc — {Qc/Q)'^ where Q denotes 
the total cms energy. 

A description of hadron spectra and hadron 
multiplicity has been previously obtained from 
the modified leading log approximation (MLLA) 
of QCD assuming a close similarity of parton and 
hadron spectra with a small transverse momen- 
tum cut-off Qo for the parton cascade (Local Par- 
ton Hadron Duality - LPHD In this de- 
scription the parameter Qo represents the con- 
finement scale which terminates the perturbative 
evolution. The multiplicity of hadrons is then 
found proportional (roughly twice in the fits) to 
the multiplicity of partons, also the energy de- 
pendences are slightly different (for a recent anal- 
ysis, see ^]). One might therefore expect that 
with decreasing fcj^ cut-off Qc the perturbative 
regime ends before the hadronic regime is reached 
at Qc Qo- 

In our analysis we search for the existence of 
such an additional scale. Furthermore we inves- 
tigate whether we can see from the experimental 
data some characteristic features of the perturba- 
tive analysis, in particular, the strong variation of 
the coupling as{k±) when approaching the QCD 
scale A. 

2. PERTURBATIVE ANALYSIS 

The inclusive and exclusive characteristics of 
the multiparticle final states can be conveniently 
treated using the generating functional. It fulfils 
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Figure 1. Comparison of different approximations to the master equation for jet evolution |g| as a function 
of the jet energy variable Y = ln(K/Qo) ('« = 2i?sin(8/2), O = tt/2) : (a) Multiplicity in quark jets; 
our exact numerical solution 0, the MLLA results from the coupled eqn. for q&g jets the solution 
for g jets only, multiplied by 4/9 and the limiting spectrum with A = 0, multiplied by 4/9 (b) 
Multiplicity ratio for exact solution ffl compared with asymptotic value 9/4 [||, the MLLA results using 
the normalization at threshold (DO ||| and CDFW and asymptotic expansions (GM and DN ^) 
for A = 0.02. 



a coupled system of evolution equations ("mas- 
ter equation") in the virtuality k = 2E sm{Q/2) 
with Q the opening angle and E the energy of the 
quark or gluon jet |^]. It takes into account the 
angular ordering inside the parton cascade as well 
as the running of the coupling with the transverse 
momentum of the emitted parton. For initializa- 
tion of the jet evolution one takes the condition 
that at the threshold of the process there is only 
one parton in the jet which carries the full mo- 
mentum. After appropriate functional differenti- 
ation one obtains the coupled integro-differential 
equation for quark and gluon jet multiplicities J\fq 
andAAg §. 

The master equation yields the complete 
MLLA asymptotics, i.e. the correct leading and 
next-to-leading terms in an expansion of y/cx^. 
Furthermore it generates a perturbative expan- 
sion which allows for the fulfilment of the initial 
conditions at threshold. The master equation can 
be simplified by making high energy approxima- 
tions already in the integral kernel ( "MLLA evo- 



lution equation" l^). For this simplified equation 
analytic expressions can be obtained for the mul- 
tiplicities in quark and gluon jets Because 
of the high energy approximations the behaviour 
near threshold has some unphysical features (mul- 
tiplicities smaller than one). To avoid this diffi- 
culty we have therefore solved the master equa- 
tion numerically for given parameters A and Qc 
(or Qo)- 

In Fig. 1 we compare the results for the mul- 
tiplicities in quark and gluon jets for different 
approximations to the master equation choos- 
ing a small value A = 0.02 for the parameter 
A = ln{Qo/A) as obtained below from a fit to 
the hadron multiplicities. One can see in (a) that 
the quark jet multiplicities for the exact solution 
and the MLLA at the same A behave very similar 
already above the inelastic threshold; the result 
for the "limiting spectrum" at A = [|| which has 
the limit A/'q ^ at threshold F — > is smaller 
in absolute normalization by roughly a factor 2. 
It should be noted however that the close agree- 
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ment within 20% of the upper three curves is spe- 
cial for A = 0.02 and disappears with increasing 
A. In Fig. lb we compare our exact solution of 
the master equation for the ratio of multiplicities 
in quark and gluon jets, again at A = 0.02, with 
several approximations using and not using the 
initial conditions at threshold. One can see that 
at LEP energies {Y sa 5) the exact solution is 
smaller than all quoted approximations. 

A qualitative difference between the exact so- 
lution of the master equation and the analytic 
MLLA results is found in the behaviour for A ^ 0: 
whereas the multiplicity approaches a finite value 
in the latter case, it diverges for the exact solution 
as it does in case of the double log approximation 
(DLA) with H 

Af ^ Ko{A^,)h{A^ln{K/A)) + ... (1) 

which at high energies and small Ac behaves like 
In Ac exp(A-\/ln(K/A)) and exhibits a logarithmic 
singularity at Ac = 0. 

The behaviour near threshold can be improved 
by taking into account the full perturbative re- 
sult at fixed order (see, for example, ref. Q). 
Such results depend on the particular process. 
For 6+ e~ annihilation the 0{as) corrected solu- 
tion has been calculated ^ as 

where Af^~'^'^*{yc) is the numerically integrated 
3 jet cross section (for running as{k±)) and 
2J\fq^\yc) is the corresponding quantity in 0{as) 
from the master equation. 

3. COMPARISON WITH e+e^ DATA 

In Fig. 2 we show the data of the average jet 
multiplicity at Q = 91 GeV as a function of the 
resolution parameter j/c obtained with the k± al- 
gorithm. The theoretical predictions from (|^) for 
the jet data are given in absolute normalization in 
terms of the single parameter A. Also shown are 
the data on hadron multiplicities in the energy 
range Q — 1.6. . . 91 GeV taken as Afaii — ^Nch 
as a function of the same scale parameter, now 
calculated as t/c = Qq/Q^, where the parameter 
Qo corresponds to the parton cutoff charac- 
terizing a hadronic scale according to the LPHD 
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Figure 2. Experimental data of the average jet 
multiplicities at LEP-1 as a function of the reso- 
lution parameter j/c! also shown are the average 
hadron multiplicities (assuming M = '^Mch) at 
different cms energies using Qc = 0.507 GeV. The 
solid (dashed) line shows the prediction in abso- 
lute normalization for the hadron (jet) multiplic- 
ities obtained by using eq. (H) with parameters 
from eq. (^ (figure from ref. |0]). 



picture and is obtained from a fit to the data. 
Another adjustable parameter here is the over- 
all normalization Kaii which relates the parton 
and hadron multiplicities in Mail = KaiiAf^orr ■ 
A good description of the data is obtained with 
parameters 

Kail = 1, A = 500±50MeV, A = 0.015±0.005 (3) 

(A = ln((5o/A)) which correspond to the curves 
in Fig. 2. The fits describe the data within 5% 
or within the errors. 

A remarkable result of our analysis with im- 
proved accuracy is the common normalization 
Kail — 1 which is possible without difficulty; the 
normalization parameter Kaii is correlated with 
Qo and can be varied within about 30%. This 
differs from the earlier results based on the limit- 
ing spectrum (approximate solution with A = 0) 
which lead to the larger value Kaii ~ 2 (see, 
for example, refs. !§,§). This difference can 
be traced back to the different normalization at 
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threshold, see Fig. la. The solution with Kaii = 
1 is natural as it provides the correct boundary 
conditions A/" = 2 at threshold for both hadrons 
and jets. Then a continuous connection between 
jet and hadron multiplicities results: both are de- 
scribed by the same equation (H) and one is ob- 
tained from the other at fixed yc by changing the 
absolute scale Qc down to Qo- There is no ev- 
idence for another nonperturbative scale above 
Qo- In this description a single hadron corre- 
sponds to a single parton of low virtuality Qo- 

An important role in this interpretation of the 
data is played by the running of as- First we 
note that in a model with fixed the multiplic- 
ity would depend only on the ratio of the available 
scales through the variable jjc but not on the ab- 
solute scales Q or Qc separately; such a model 
would predict at high energy a power like depen- 
dence on Q/Qc, i-e. a straight line in Fig. 2 for 
both the hadron and jet multiplicities between 
the two curves shown. 

On the other hand, in case of running as the 
additional scale kj_/ A appears which, for a given 
He, leads to an enhancement of multiplicities at 
small absolute scales Q or Qc and hence to the 
different curves in Fig. 2. The effect of the run- 
ning as is most pronounced at large yc- Here 
the quantity A/" — 2 is of 0{as) and therefore re- 
flects the typical coupling. From the separation 
of the jet and hadron data in Fig. 2 (see also data 
on AA — 2 in one finds that around yc ~ 0.1 
the typical couplings should differ by more than 
an order of magnitude. As the coupling to pro- 
duce jets at LEP energies is around as ~ 0.1 the 
coupling to produce particles at low energies is 
therefore 



> 1 at Q - 1.6 GeV. 



(4) 



The fact that both data sets are rather well de- 
scribed over the full kinematic region where data 
are available confirms that kj_ is a good choice 
for the scale of as and also supports the LPHD 
concept to describe hadrons by the k± cut-off Qo- 
There is a certain ambiguity in the definition of 
k± which leads to different values for A. The 
value quoted above in (^) refers to the Durham- 
fc_L definition whereas the usual k± would yield a 
smaller value around 300 MeV [pJ. 



Another characteristic sign of the running cou- 
pling is the strong rise of the jet rate towards 
small j/c- Only for very high resolution, if Qc 
is lowered from 900 MeV {yc 10"") to the fi- 
nal Qo ^ 500 MeV, about three quarters of the 
final multiplicity are produced, three times as 
many jets as in the large complementary kine- 
matic region with j/c > lO"*^. The steep rise to- 
wards small j/c in Fig. 2 refiects the Landau sin- 
gularity in the coupling at j/a = A^/Q^, which 
is however screened by the hadronization scale 
2/0 = Qo/Q^ ^ 2/A- This behaviour is qualita- 
tively described by the high energy DLA result 
(|l|) which shows the singularity at A = 0. It will 
be interesting to study this behaviour experimen- 
tally, but some additional assumptions on mass 
effects are necessary in this region ||l[ . 

The prediction for the hadron multiplicity in 
gluon jets (Fig. lb) is also well met by experi- 
mental data at LEP energies [Q; there are also 
predictions for the dependence of the subjet mul- 
tiplicity in gluon jets on the scale parameter j/c- 
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